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ABSTRACT 

Wolf-Rayet stars represent one of the final stages of massive stellar evolution. 
Relatively little is known about this short-lived phase and we currently lack re- 
liable mass, distance, and binarity determinations for a representative sample. 
Here we report the first visual orbit for WR 140(=HD193793), a WC7+05 bi- 
nary system known for its periodic dust production episodes triggered by intense 
colliding winds near periastron passage. The IOTA and CHARA interferome- 
ters resolved the pair of stars in each year from 2003-2009, covering most of the 
highly-eccentric, 7.9 year orbit. Combining our results with the recent improved 
double-line spectroscopic orbit of Fahed et al. (2011), we find the WR 140 sys- 
tem is located at a distance of 1.67 ± 0.03 kpc, composed of a WR star with 
M W r = 14.9 ± 0.5 M and an O star with M Q = 35.9 ± 1.3 M . Our precision 
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orbit yields key parameters with uncertainties ~ x6 smaller than previous work 
and paves the way for detailed modeling of the system. Our newly measured 
flux ratios at the near-infrared H and Ks bands allow an SED decomposition and 
analysis of the component evolutionary states. 

Subject headings: stars: binaries: visual, stars: individual (WR 140, HD 193793), 
stars: Wolf-Rayet, techniques: interferometric, infrared: stars 



Introduction 



Wolf-Rayet (WR) stars are evolved massive stars characterized by intense mass-loss 
through radiation-driven winds. These hot, emission-line stars are especially helium-rich, 
having lost most of their hydrogen envelope through winds or interaction with a companion. 
The range of progenitors that become WR stars is not well understood and establishing a 
massive star evolution sequence re presents one of the most serious challenges for modern 
stellar theory (see recent review by ICrowtherl 120071 ) . 



Concrete mass and distance determinations are crucial to making further progress, but 
this is difficult due to the large distances from Earth at which most WR stars lie. According to 
van der Huchtl (120011 ). there are only 19 WR stars with mass estimates based on spectroscopic 
orbits and the vast majority of these are for short-period systems with periods between 1-100 
days. Longer period binaries are less likely to have had interactions between the components 
but are difficult to characterize due to the lower orbital speeds. 

The subject of this Letter i s WR 140 (=HD 193793), a WR binary system with a 
7.9 year period. IWilliams et al.l (119871) first noticed mysterious episodic infrared variabil- 
ity and follow-up observations ( IWilliams et al.l Il990l : iMoffat et al.l 119871 ) established the 
cause to be dust creation near periastron of a highly eccentric orbit, likely catalyzed in the 
colliding-wind interface between the WR star and O-star winds ( IUsovlll99ll ). Since this time, 
WR 140 (WC7+Q5) has been subject to many mon itoring; campaigns, inclu ding infrared 



(IWil 



iams et al.l 120091 ). radio (IWhite fc Becker! Il995l ). and in radial velocity (IFahed et al. 



201 ih . Recently, progress towards a proper visual orbit was made by the single-epoch detec- 



tion of the binary using the IOTA interferometer (jMonnier et al. 



200 4J) and through rep eated 



imaging of the rotating colliding-wind region using the VLB A (iDougherty et al.l 120051 ). De- 
spite the wealth of data, mass estimates have suffered from large errors (~20%) due to lack 
of a high-quality visual orbit to go along with precise spectroscopic data. 



In this paper, we report seven epochs of binary observations at the near-infrared H and 
Ks bands with the IOTA and CHARA interferometers. Our data span seven years allowing 
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us to construct the first complete visual orbit for WR 140. We combine this with recent 
spectroscopic work to determine precise masses and orbital parallax. 



Observations 



2.1. IOTA 



Observations of WR 140 in 2003, 200 4, and 2005 were o btained with the IOTA (Infrared- 
Optical Telescope Array) interferometer ( ITraub et al.ll2003f ). IOTA was located on Mt. Hop- 
kins (Arizona) and consisted of three 0.45-m telescopes that were movable among 17 stations 
along 2 orthogonal linear arms (telescopes A & C could move along the 35-m northeastern 
arm, while telescope B could move along the 15- m southeastern arm). By observing a target 
in many different array configurations, IOTA could synthesize an aperture 35m x 15m (corre- 
sponding to an angular resolution of ^5x 12 milliarcseconds at 1.65,um). The observations 
of WR 140 from 2003 were first reported in lMonnier et al.l (12004] ) . 



All of the IOTA observations included three simultaneous baselines using the broad 
band H filter and the light beams from the th ree telescopes were interfered using the single- 
mode IONIC3 combiner (IBerger et al.l 120031 ). Basic data reduction procedure s were the 



same as described in several previous IOTA papers fe. g.. iMonnier et al.l 120061 ) and final 



error estimation followed the study of the A Vir binary by lZhao et al.l (120071 ); we applied 2% 
relative and a AV 2 = 0.02 additive systematic errors to all our measured visibility amplitude 
V 2 data. We found our closure phases were partially corrupted by bandwidth smearing effects 
and we only used V 2 dat a for our orbit fit ting (see detailed description and simulation of 
these effects explored by IZhao et al.l 120071 ). Table [1] contains a detailed log of the IOTA 
observations. 



2.2. CHARA 



Observations of WR 140 in the years 2006, 2007, 2008, and 2009 were obtained with the 
Center for High Angular Resolution Astronomy (CHARA) Array. CHARA is located on Mt. 
Wilson (California) and consists of six fixed 1-m aperture telescopes with baselines ranging 
from 30 to 330m. CHARA is currently the longest baseline optical interferometer in the 
world and can reach angular resolutions of ~0.5 milli-arcseconds (mas) in the near-infrared. 
Th e CHARA facility and t he t wo-beam CLASSIC combiner used in this work are described 



by lten Brummelaar et al.l (120051 ) . 
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Most of the CHARA observations were carried out in the Ks band although a few data 
points were taken in the H band. The observing dates, wavelengths, and baselines can be 
found in Tabled) We reduced the data using an IDL-ba sed suite of routines writt en by one 
of us (JDM) and these have been previously described in iTannirkulam et al.l (120081 ). For this 
paper, we adopted a 10% relative and a AV 2 = 0.02 additive error to account for systematic 
calibration error^]. 

We have reported all the calibrators and adopted sizes in Table [TJ Uniform Disk 
(UD) diameters of interferometer calibrators were generally estimated using getCal, an 
SED-fitting routine maintained and distri buted by the NASA Exoplanet Science Institute 
(http:/ /nexsci.caltech.edul), or SearchCal ( IBonneau et al.l 120061 ). a surface brightness algo- 



rithm maintained and distributed by the J. M. Mariotti Center (JMMC, http://jmmc.fr). 
We note that getCal does not produce as accurate diameter estimates as searchCalm general, 
but getCal can be employed on a wider variety of calibrator spectral types and was more 
useful in this current paper. 



The final IOTA and CHARA calibrated V 2 data were saved in OI-FITS format ( jPauls et al 



20051 ) and all files are available upon request. 



3. Analysis 

We have interpreted our V 2 data using a binary model. Here we describe the components 
of the model and our fitting results. 



3.1. Description of binary model 

Our model for WR 140 consists of two stars, each modelled as a uniform disk. The 
size of the O-star can be estimated based on the effective temperature and infrared flux 
estimate, wh ile the WR-star size is affected by the optically-thick wind (we followed similar 
procedure as lMillour et al.ll2007l . to account for wind opacity). We found that both stars are 
much smaller than the resolution of CHARA and so our final results are not sensitive to our 
adopted UD size of 0.05 mas for the WR star and 0.07 mas for the O-star. 



The WR/O-star flux ratios at H band and Ks band were separately fitted but held 



*The calibration errors were lower for IOTA-IONIC3 because of the use of single- mode fibers while 
CHARA-CLASSIC is a "free-space combiner" with only limited spatial filtering. 
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const ant across all epochs. The near-infrared flux monitoring by iTaranova fc Shenavrin 



( 120 111 ) showed that 52% of the K-band emission during the 2009 epoch was from the outburst 
dust shell and was assumed to be over-resolved by CHARA. No other epochs were affected 
by dust emission. For a given set of orbital elements (a, e, i, cu, Q, P, T ), we can predict the 
separation and position angle at the time of each observation. The apparent brightness ratio 
can be affected by the finite bandwidth of the observations ( "bandwidth-smearing" ) and we 
have account ed for the square bandpasses of the H and Ks filter using the basic procedure 



described in IZhao et al.l ( 120071 ). The bandwidth-smearing correction was insignificant for 
the short baseline IOTA data, but did affect the longer baseline CHARA data at the 5-10% 
level. 

For visualization purposes, we also collected each year's data and fitted for relative 
positions. The best fitting locations of the O-star relative to the WR star for each observation 
year has been included in Figured] and allows comparison with our final orbit fitting results. 
The size and shape of each epoch's allowed region (an error ellipse containing the 68% 
confidence region) vary significantly from year to year due to differences in the quantity and 
quality of the V 2 data. We also included the parameters of these best-fit locations in Table HJ 

To validate our choice of this simple binary model, we also carried out model-independent 
image reconstructions for the 2003-2005 IOTA/IONIC3 data. We confirmed that the system 
is dominated by two point sources. We were interested to see if there could be any sign of 
the colliding wind zone between the stars, but no extra emission was seen in this region (at 
the level of a few percent of the peak emission). The brightest emission not coming from 
the two stars showed up to the North-East of the system in some epochs at the 2-5% level 
- likely an artifact from residual miscalibration. 



3.2. Orbit fitting procedure 



To arrive at our final orbital solutions, we fit directly to the V 2 values. We carried out 
two different fitting exercises that differed in how we incorporated spectroscopic data. 

First we wanted to carry out an or bital fit as independent as possible from the recent 
spectroscopic orbit of lFahed et al.l ( 120111 ) . This allows us to independently confirm the crucial 
orbital elements e and u, although we adopted their values for the period (P) and time of 
periastron (T ) in this procedure. The best-fitting orbital elements (reduced x 2 — 0.52) are 
compiled in t he right column of Table [3j Error bars were estimated using 1000 bootstrap 
resamplings (jEfron &: Tibshiranil Il993l ) of the data (group ed by night). This is the same 



procedure recently applied for the visual orbit of a Oph (IHinkley et al.l 1201 ll ) . In order 
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to capture the uncertainties in the iFahed et al.l (1201 1 ) estimates for P and To, we did not 
strictly fix th ese quantitie s durin g the bootstrap fits but rather used Monte Carlo sampling 
based on the IFahed et al.l (120111 ) uncertainty estimates. 



The IOTA+CHA RA visual orbit com pares favorably with the spectroscopically-determined 
orbital elements from IFahed et al.l (120111 ): 



e: 0.8962±0.0014 (Fahed's orbit "This paper+M03") compared to eccentricity 0.90ll;^ 
(this work alone). These values are compatible and confirm the high orbit eccentricity. 

w: 44.6°±1.1° (Fahed's orbit "This orbit+M03") compared to 48.2°±1.3°(this work 
alone). This is slightly discr e pant , although note that one of the orbital solutions 
presented in the IFahed et al.l ( 1201 ll ) report u = 47.5° (Fahed's orbit "This paper"; 
middle solution in Table 2). 



Much tighter co nstraints on t he or bital elements can be attained by using the spec- 
troscopic values from IFahed et al.l ( 1201 ll ) as a prior in the visual orbit fit (the orbit labeled 
"This paper+M03"). The m ost importan t effec t of this is to constrain the eccentricity, which 
is better determined by the IFahed et al.l ( 1201 ll ) dataset that sampled the periastron period 
very densely. Our visual observations missed the fast-changing orbital motion in early 2009 
and thus can not be expected to optimally constrain eccentricity. 

The first column of Table [3] contains our best-fitting solution from the joint analysis 
(reduced \ 2 — 0.51). We used the same bootstrap procedure to determine the error bars 
on each parameter. To visualize the range of orbits allowed by our solutions, we plotted all 
1000 bootstrap orbits in Figure HJ As expected, the joint solution is better constraining, 
especially near periastron. 

The last step in our analysis was to take the new orbital elements and use the spectro- 
scopic K values from IFahed et al.l (120111 ) to calculate masses and orbital parallax. Strictly 
speaking, since the K values depend slightly on the orbital elements themselves, we re-fitted 
the 7 and K values for our orbital solutions using the radial velocity data in IFahed et al. 
(120111 ). This refinement is slightly more accurate than simply adopting the spectroscopic 
a sin i along with astrometry to derive inclination and distance. The final results for masses 
and distance are also included in Table [3j 

Another way to view the remarkable quality of the fit is presented in Figure [2J Here we 
show the observed V 2 as a function of the projected separation of the two components (joint 
orbital solution). This figure also shows the effect of bandwidth smearing for large projected 
separations. 
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Discussion and Future Work 



The WR and O-star masses found here are similar to previous estimates (e.g. Ivan der Hucht 

20011 ) but with about 6x smaller error bars (~3-4%). As has been seen in other WC bi- 



nary systems, we find the WR star mass is less than 20 M and less than half the current 
mass of the O-star companion. These two qualit ies fit the interesting trend seen in the (six) 
WC stars with mass estimates (jCrowtherl 120071 ) . Our accurate distance will help to place 
these two stars in the Hertzsprung-Russel diagram for a stringent test of stellar evolution 
models that include mass-loss. This is only the second galactic WR star to have a distance 
measured through (orbital) parallax (the other is 7 Vel; iMillour et al.l 120071 ) . Note that the 
periastron distance of the orbit is 1.53 AU, which means the O-star companion is too distant 
to have affected the stellar evolution of the Wolf-Rayet s tar - unless the WR star progenitor 
experienced a red supergiant stage (see case of WR 104; iTuthill et al.ll2008l ) or the orbit has 
drastically changed due to mass-loss. 

. Fluxr 



The interferometry data allow us to measure the flux ratios of the two stars 



"(WR) ■ 



Flux 



(O) 



in the H- and Ks-bands for the first time. Our H-band IOTA data employed closure phases, 
allowing us to identify the Northwest component to be brightest in 2003. We have had to 
assume that this component is also brighter at K band since we lack closure phase data to 
break the 180° degeneracy of single-baseline data. 

Based on our flux ratios, the IR-bright component has a significantly redder color than 
the IR-faint component, consistent with the IR-bright compon ent being the WR star. The 
shape of the non-thermal radio emission (jDougherty et al.ll2005f ) also identifies the Northwest 
(IR-bright) component in 2003 to be the WR star that is expected to possess the the higher- 
momentum wind. 

Armed with this knowledge, we can decompose the combined spectral energy distribu- 
tion (SED) into their component SEDs. For the combined system, we adopt (non-dusty) 



BVJHK magnitudes = (7.28, 6.89. 5.71. 5.35. 5. 02) flReedll2003l: iTaranova fc Shenavrinll201lh 



Morris et al.l fll993h 



1990 



and assume the 05III star has colors given by iMartins fc Plezl (120061 ). The last ingredient 
w e need i s the l east certain: the interstellar reddening. We will adop t the reddening law 
of iMathid (ll990l ) with Ry=3A and with Ay spanning A y = 2.95 from 
using the 2175A feature and A v = 2.06 from SED colors ( Conti Sz Vacca 
see, the large uncertainty in Ay leads to large errors in our system luminosity despite our 
new well-constrained distance. We can now fix the H band ratio to 1.37 from our IOTA data 
and solve for the remaining flux ratios (modulo the Ay uncertainty). For Ay = 2.95, the flux 



- as you will 



'''Here, we have used the relation that A v = 1.1 Ay (e.g., Smith 19681 ) 
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ratios at BVJHK become (2.46, 1.64, 1.05, 1.37, 1.93); for Ay = 2.06 the flux ratios become 
(0.37,0.35,0.88,1.37,2.09). The observed flux ratio at Ks band (1.94±0.06) is within the 
range seen here and is in good agreement with the high Ay ca se. Depending on which lines 
are used and which templates are adopted, iFahed et al.l (120111 ) argue that the O-star might 
range from ~0.5 to ~3x brighter than the WR star in the optical band continuum - a wide 
range compatible with our derived flux ratios. Our prediction for V band flux rato varies 
by a factor of 5 (!) because of Ay uncertainties and thus an interferometric measurement 
here could be exploited to strongly constrain the true reddening, a crucial component to 
luminosity determination which wei now explore. 

Using our new distance (d = 1.67 kpc) and the range of flux ratios, we find My for 
the O-star ranges from -6.11 and -5.94 and the My for the WR star ranges from -6.6 and 
-4.8. This would cla ssify the 05 star as intermediate between giant and supergiant (cf. 
Martins fc Plezl 120061 ). The WR star luminosity class is hardly constrained but the high Ay 
case yields unrealisticly hi gh luminosities as jud ged by other WC7 stars with known distances 
(via cluster membership; Ivan der Huchtl 120011 ). In an attempt to recon cile the diffe r ences , 
we scaled the fluxes for a 70kK CMFGEN model for the WC7 star from lSmith et al.l J2OO2I ) 
to m atch the measured flu x ratio of 1.37 in the H Band using a model SED for the 05 star 
from iMartins et al.l ( 120051 ). The WR/05 flux ratios in V and v were found to be 1.05 and 
0.73 respectively, that in V being affected by strong emission lines. The absolute magnitudes 
derived for the 05 star, My = -6.37 and -5.65 for high and low reddening extremes, are near 
or above those for supergiants, whereas those for the WC star, M v = -6.7 and -5.9 are both 
anomalously high. Resolution of the uncertain interstellar reddening of WR 140 is beyond 
the scope of the present paper but is urgently required to exploit the determination of its 
parallax. 

We expect our success here will motivate future observations of galactic Wolf-Rayet 
stars using today's interferometers from visible to NIR wavelengths. The mult i- wavelength 
flux ratios can help yield crucially-needed new estimates of interstellar reddening and the 
possibility to combine spectroscopic orbits with new interferometric visual orbits will al- 
low accurate distance, mass, and luminosity measurements for a substantially larger set of 
galactic WR stars. Such a dataset is important to test the current massive stellar evolu- 
tion paradigm that tells us how main-sequence O-stars move through the various stages of 
Red Supergiant, Wolf-Rayet (WN & WC) and Luminous Blue Variable before ultimately 
becoming a Supernova. 



We have appreciated discussions with Tony Moffat, Peter Tuthill, Debra Wallace, Bill 
Danchi, Sean Dougherty, and Remi Fahed during the (long) course of this work. We thank 
SAO, U. Mass, NSF AST-0138303, NSF AST-0352723, and NASA NNG05G1180G for sup- 



-9 - 



porting IOTA development and operations. We also acknowledge funding from GSU, the 
Keck Foundation, and NSF AST-0908253 for the CHARA Array. IONIC-3 was developed 
by LAOG (now IPAG) and LETI in the context of the IONIC collaboration (LAOG, IMEP, 
LETI), funded by the CNRS and ONES (France). Lastly we thank NSF AST-0807577 
for support of University of Michigan researchers in this work. EP received funding from 
a Michelson Postdoctoral Fellowship and a Scottish Universities Physics Alliance (SUPA) 
advanced fellowship. PMW is grateful to the Institute for Astronomy for hospitality and 
continued access to the facilities of the Royal Observatory Edinburgh. This research has 
made use of the SIMBAD database, operated at CDS, Strasbourg, France, and NASA's 
Astrophysics Data System (ADS) Bibliographic Services. 

Facility: IOTA (IONIC3), CHARA (CLASSIC) 



REFERENCES 

Berger, J., Haguenauer, P., Kern, P. Y., Rousselet-Perraut, K., Malbet, F., Gluck, S., Lagny, 
L., Schanen-Duport, I., Laurent, E., Delboulbe, A., Tatulli, E., Traub, W. A., Car- 
leton, N., Millan-Gabet, R., Monnier, J. D., & Pedretti, E. 2003, in Interferometry for 
Optical Astronomy II. Edited by Wesley A. Traub . Proceedings of the SPIE, Volume 
4838, pp. 1099-1106 (2003)., 1099-1106 

Bonneau, D., Clausse, J.-M., Delfosse, X., Mourard, D., Cetre, S., Chelli, A., Cruzalebes, 
P., Duvert, G., & Zins, G. 2006, A&A, 456, 789 

Conti, P. S. & Vacca, W. D. 1990, AJ, 100, 431 

Crowther, P. A. 2007, ARA&A, 45, 177 

Dougherty, S. M., Beasley, A. J., Claussen, M. J., Zauderer, B. A., & Bolingbroke, N. J. 
2005, ApJ, 623, 447 

Efron, B. & Tibshirani, R. J. 1993, An introduction to the bootstrap (New York: Chapman 
and Hall) 

Fahed, R., Moffat, A. F. J., Zorec, J., Eversberg, T., Chene, A.-N., Alves, F., Arnold, W., 
Bergmann, T., Carreira, L. F. G., Dias, F. M., Fernando, A., Gallego, J. S., Hunger, 
T., Knapen, J. H., Leadbeater, R., Morel, T., Rauw, G., Reinecke, N., Ribeiro, J., 
Romeo, N., dos Santos, E. M., Schanne, L., Stahl, O., Stober, B., Stober, B., Viegas, 
N. G. C, Vollmann, K., Corcoran, M. F., Dougherty, S. M., Pittard, J. M., Pollock, 
A. M. T., & Williams, P. M. 2011, MNRAS, in press 



- 10 - 



Hinkley, S., Monnier, J. D., Oppenheimer, B. R., Roberts, Jr., L. C, Ireland, M., Zimmer- 
man, N., Brenner, D., Parry, I. R., Martinache, F., Lai, O., Soummer, R., Sivara- 
makrishnan, A., Beichman, C, Hillenbrand, L., Zhao, M., Lloyd, J. P., Bernat, D., 
Vasisht, G., Crepp, J. R., Pueyo, L., Shao, M., Perrin, M. D., King, D. L., Bouchez, 
A., Roberts, J. E., Dekany, R., & Burruss, R. 2011, ApJ, 726, 104 

Martins, F. & Plez, B. 2006, A&A, 457, 637 

Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 

Mathis, J. S. 1990, ARA&A, 28, 37 

Millour, F., Petrov, R. G., Chesneau, O., Bonneau, D., Dessart, L., Bechet, C, Tallon- 
Bosc, I., Tallon, M., Thiebaut, E., Vakili, F., Malbet, F., Mourard, D., Antonelli, 
P., Beckmann, U., Bresson, Y., Chelli, A., Dugue, M., Duvert, G., Gennari, S., 
Gliick, L., Kern, P., Lagarde, S., Le Coarer, E., Lisi, F., Perraut, K., Puget, P., 
Rantakyro, F., Robbe-Dubois, S., Roussel, A., Tatulli, E., Weigelt, G., Zins, G., 
Accardo, M., Acke, B., Agabi, K., Altariba, E., Arezki, B., Aristidi, E., Baffa, C, 
Behrend, J., Blocker, T., Bonhomme, S., Busoni, S., Cassaing, F., Clausse, J.-M., 
Colin, J., Connot, C, Delboulbe, A., Domiciano de Souza, A., Driebe, T., Feautrier, 
P., Ferruzzi, D., Forveille, T., Fossat, E., Foy, R., Fraix-Burnet, D., Gallardo, A., 
Giani, E., Gil, C, Glentzlin, A., Heiden, M., Heininger, M., Hernandez Utrera, O., 
Hofmann, K.-H., Kamm, D., Kiekebusch, M., Kraus, S., Le Contel, D., Le Contel, 
J.-M., Lesourd, T., Lopez, B., Lopez, M., Magnard, Y., Marconi, A., Mars, G., 
Martinot-Lagarde, G., Mathias, P., Mege, P., Monin, J.-L., Mouillet, D., Nussbaum, 
E., Ohnaka, K., Pacheco, J., Perrier, C, Rabbia, Y., Rebattu, S., Reynaud, F., 
Richichi, A., Robini, A., Sacchettini, M., Schertl, D., Scholler, M., Solscheid, W., 
Spang, A., Stee, P., Stefanini, P., Tasso, D., Testi, L., von der Liihe, O., Valtier, 
J.-C, Vannier, M., & Ventura, N. 2007, A&A, 464, 107 

Moffat, A. F. J., Lamontagne, R., Williams, P. M., Horn, J., & Seggewiss, W. 1987, ApJ, 
312, 807 

Monnier, J. D., Berger, J.-P., Millan-Gabet, R., Traub, W. A., Schloerb, F. P., Pedretti, E., 
Benisty, M., Carleton, N. P., Haguenauer, P., Kern, P., Labeye, P., Lacasse, M. G., 
Malbet, F., Perraut, K., Pearlman, M., & Zhao, M. 2006, ApJ, 647, 444 

Monnier, J. D., Traub, W. A., Schloerb, F. P., Millan-Gabet, R., Berger, J.-P., Pedretti, 
E., Carleton, N. P., Kraus, S., Lacasse, M. G., Brewer, M., Ragland, S., Ahearn, A., 
Coldwell, C, Haguenauer, P., Kern, P., Labeye, P., Lagny, L., Malbet, F., Malin, D., 



- 11 - 



Maymounkov, P., Morel, S., Papaliolios, C, Perraut, K., Pearlman, M., Porro, I. L., 
Schanen, L, Souccar, K., Torres, G., & Wallace, G. 2004, ApJ, 602, L57 

Morris, P. W., Brownsberger, K. R., Conti, P. S., Massey, P., & Vacca, W. D. 1993, ApJ, 
412, 324 

Pauls, T. A., Young, J. S., Cotton, W. D., & Monnier, J. D. 2005, PASP, 117, 1255 
Reed, B. C. 2003, AJ, 125, 2531 
Smith, L. F. 1968, MNRAS, 140, 409 

Smith, L. J., Norris, R. P. F., & Crowther, P. A. 2002, MNRAS, 337, 1309 

Tannirkulam, A., Monnier, J. D., Harries, T. J., Millan-Gabet, R., Zhu, Z., Pedretti, E., Ire- 
land, M., Tuthill, P., ten Brummelaar, T., McAlister, H., Farrington, C, Goldfinger, 
P. J., Sturmann, J., Sturmann, L., & Turner, N. 2008, ApJ, 689, 513 

Taranova, O. G. & Shenavrin, V. I. 2011, Astronomy Letters, 37, 30 

ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., Bagnuolo, Jr., W. G., Turner, 
N. H., Sturmann, L., Sturmann, J., Berger, D. H., Ogden, C. E., Cadman, R., 
Hartkopf, W. L, Hopper, C. H., & Shure, M. A. 2005, ApJ, 628, 453 

Traub, W. A., Ahearn, A., Carleton, N. P., Berger, J., Brewer, M. K., Hofmann, K., Kern, 
P. Y., Lacasse, M. G., Malbet, F., Millan-Gabet, R., Monnier, J. D., Ohnaka, K., 
Pedretti, E., Ragland, S., Schloerb, F. P., Souccar, K., & Weigelt, G. 2003, in Inter- 
ferometry for Optical Astronomy II. Edited by Wesley A. Traub. Proceedings of the 
SPIE, Volume 4838, pp. 45-52 (2003)., 45-52 

Tuthill, P. G., Monnier, J. D., Lawrance, N., Danchi, W. C, Owocki, S. P., & Gayley, K. G. 
2008, ApJ, 675, 698 

Usov, V. V. 1991, MNRAS, 252, 49 

van der Hucht, K. A. 2001, New A Rev., 45, 135 

White, R. L. & Becker, R. H. 1995, ApJ, 451, 352 

Williams, P. M., Marchenko, S. V., Marston, A. P., Moffat, A. F. J., Varricatt, W. P., 
Dougherty, S. M., Kidger, M. R., Morbidelli, L., k Tapia, M. 2009, MNRAS, 395, 
1749 



- 12 - 



Williams, P. M., van der Hucht, K. A., Florkowski, D. R., Pollock, A. M. T., & Wamsteker, 
W. M. 1987, in IAU Symposium, Vol. 122, Circumstellar Matter, ed. I. Appenzeller 
& C. Jordan, 453-454 

Williams, P. M., van der Hucht, K. A., Pollock, A. M. T., Florkowski, D. R., van der Woerd, 
H., & Wamsteker, W. M. 1990, MNRAS, 243, 662 

Zhao, M., Monnier, J. D., Torres, G., Boden, A. F., Claret, A., Millan-Gabet, R., Pedretti, 
E., Berger, J. -P., Traub, W. A., Schloerb, F. P., Carleton, N. P., Kern, P., Lacasse, 
M. G., Malbet, F., & Perraut, K. 2007, ApJ, 659, 626 



This preprint was prepared with the A AS IATgX macros v5.2. 



13 



Table 1. Observing log for WR140. All calibrated OI-FITS data available upon request. 



Orbital 


Date 


Phase a 


CUT) 


2. 


.296 


2003Junl7 


2 


.406 


2004Apr30 


2 


.407 


2004May01 


2 


.417 


2004May30 


2 


.417 


2004Jun01 


2 


.418 


2004Jun04 


2 


.419 


2004Jun05 


2 


.419 


2004Jun06 


2 


.547 


2005Junll 


2 


.547 


2005Junl3 


2 


.548 


2005Junl4 


2 


.548 


2005Junl5 


2 


.548 


2005Junl6 


2 


.549 


2005Junl7 


2 


.549 


2005Junl8 


2 


.698 


2006Aug22 


2. 


.799 


2007Junll 


2 


.800 


2007Junl3 


2 


.926 


2008Junl4 


2. 


.927 


2008Junl6 


2. 


.927 


2008Junl7 


2 


.928 


2008Junl8 


3 


.054 


2009Jun20 


3 


.055 


2009Jun22 


3 


.055 


2009Jun23 


3 


.056 


2009Jun25 



Interferometer 
(Configuration) 



IOTA b (A35C15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
IOTA (A35B15C10) 
CHARA C (W2-E2) 
CHARA (W2-E2) 
CHARA (W2-E2) 
CHARA (W2-E2) 
CHARA (W2-S1) 
CHARA (W1-S2) 
CHARA (W1-S2) 
CHARA (W1-S2) 
CHARA (W1-S2) 
CHARA (S2-E2) 
CHARA (S2-E2) 



Ao Bandwidth 
(H Om) 



1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


1 


.650 


0.248 


2 


.133 


0.350 


1 


.673 


0.274 


2 


.133 


0.350 


1 


.673 


0.274 


2 


.133 


0.350 


2 


.133 


0.350 


2 


.133 


0.350 


2 


.133 


0.350 


2 


.133 


0.350 


2 


.133 


0.350 


2 


.133 


0.350 



a Orbital phase assuming T = 2446156.2 (MJD), P = 2896.5 days 
iFahed et alll201ll'l . 

b IOTA calibration employed the following calibrators (all sizes were esti- 
mated using getCal): HD 192985 (0.46±0.06 mas), HD 193631 (0.31±0.28 
mas), HD 126035 (0.78±0.24 mas), HD 193664 (0.58±0.05 mas), HD 193961 
(0.24±0.06 mas) 

c CHARA calibration employed the following calibrators (all sizes es- 
timated using getCal, except HD196360): HD 192985 (0.46±0.06 mas), 
HD 193631 (0.31±0.28 mas), HD 196360 (0.61±0.05 mas), HD 192640 
(0.40±0.15 mas), HD 195194 (0.63±0.13 mas) 
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Table 2. Position of O-star with respect to WR star 



Mean Date 


Relative Position 




Error Ellipse a 


(UT) 


East (mas) 


North (mas) 


Major (mas) 


Minor (mas) 


PA Major (E of N) 


2003Junl7 


5.70 


-11.31 


0.87 


0.30 


-86 


2004May25 


6.87 


-11.33 


0.34 


0.11 


-61 


2005Junl5 


7.36 


-10.16 


0.37 


0.11 


-57 


2006Aug22 


6.97 


-7.14 


4.85 


0.26 


112 


2007Junl2 


5.73 


-4.18 


1.25 


0.10 


-28 


2008Junl7 


3.51 


-0.76 


0.11 


0.05 


51 


2009Jun23 


0.86 


-5.46 


0.70 


0.13 


-47 



a Error ellipse contains 68% confidence interval and is specified by the ±error in the two orthogonal 
directions (in milliarcseconds) specified by the position angle of the ellipse major axis (degrees East of 
North). See Figure [T] for graphical representation. 
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Table 3. Orbital parameters for WR140 (adopting WR as primary star, as secondary 

star) 



Parameter 



This work with 
Fahed et al. (2011) prior 



This work alone a 



Flux Ratio b (H band) 
Flux Ratio (Ks band) 
Semi-Major Axis (mas) 
Eccentricity 
Inclination (deg) 
u (deg) 
Q (deg) 
Period (days) 
T (MJD) 

ExVdof 



1.37±0.03 
1.94±0.06 



8.82±0.05 

0.8964±g : gSS 
119.6±0.5 
46.8±0.4 
353.6±0.4 
2896.35±0.20 
46154.8±0.8 
0.51 



8 qq+o.16 



o.90i; 



.006 
.004 

118.9toj 
48.2±1.9 

354.218:1 
(2896.5±?:§) 

+2.6\ 
3.3J 



(46155.7 
0.52 



Derived Physical Quantities 


Distance (kpc) 


1.67±0.03 


i fin+o.ii 
i ' ou -0.07 


M W r (M ) 


14.9±0.5 


13.9ii| 


M (M ) 


35.9±1.3 


33.1+IJ 



a All orbital elements were fitted to the visual orb it date presented he re, 
except for the P and To, which were adopted from iFahed et al.l ( 1201 ll ). 



b FluX ratio is WR flux density _ 
(J5 flux density 



The visual orbit above and the spectroscopic data from IFahed et al. 



(120111 ) were combined to derive the orbital parallax and relative masses. 
We emphasize that the high precision on the mass and distance requires 
b oth the visua l orbit presented here and high quality spectroscopic data 
of IFahed et all (boilh . 
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Fig. 1. — This figure shows the 3-ct band of allowed orbits based on our visual data alone (thick blue band) 
and for the more-constraining, joint visual/spectroscopic solution (yellow band). The best fit joint solution 
is shown with solid line. For each year, the data were analyzed separately and the positions of the O star 
with respect to WR star are shown here, marked by error ellipses (see Table [2]). The best-fit orbit prediction 
for each epoch is connected to each correponding error ellipse. 
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Fig. 2. — Here we plot the observed visibility data as a function of the projected binary separation of 
our best-fitting joint model. We have plotted the expected curve both for no bandwidth smearing (solid 
line) and for the appropriate level we used in this work (dashed line). The bottom panels shows the residuals 
between data and model normalized by the data errors. Note the visibility for the CHARA 2009 data was 
boosted by factor of 2.1 to account for the ext ra over-resolved emission by t he short-lived dust shell created 
at periastron (48% stellar emission, 52% dust; iTaranova fc Shenavrinll201ll) . 



